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ABSTACT 

The  in-coupling  process  for  grating-coupled  planar  optical  waveguide  sensors  is  investigated 
in  the  case  of  TE  waves.  A  simple  analytical  model  based  on  the  Rayleigh-Fourier  method  is 
applied  together  with  a  perturbational  technique  to  calculate  analytical  expressions  for  the 
guided  wave  amplitudes.  In  addition,  analytical  expressions  are  derived  for  the  position 
correction  and  width  of  the  in-coupling  resonant  peaks.  Numerical  computations  verify  the 
model  for  shallow  gratings  both  in  terms  of  peak  shape  and  position  and  provide  the 
limitations  for  the  analytical  formulas. 
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1.  Introduction 


The  application  of  grating  couplers  in  planar  optical  waveguides  was  introduced  by  Lukosz  and 
Tiefenthaler1  in  1983  as  transducer  elements  for  biological  and  chemical  sensing.  Since  the 
introduction  of  this  concept  a  vast  amount  of  scientific  effort  has  been  invested  into  the 
experimental  and  theoretical  characterization  of  the  grating  coupler  in  terms  of  optimizing  its 
sensitivity2"8. 

In  most  cases,  the  grating  coupler  consists  of  a  sinusoidal  or  square-wave-formed  surface 
corrugation  embedded  in  the  waveguiding  film,  implying  that  the  actual  thickness  of  the  film 
varies  periodically  in  the  grating  region.  This  thickness  modification  alters  the  mode  properties  of 
the  waveguide  as  compared  to  those  of  a  non-corrugated  waveguide.  In  the  literature,  this 
influence  is  often  neglected,  as  the  grating  amplitude  is  considered  small  in  comparison  with  the 
film  thickness.9"14  However,  Kunz  et  al.15  exposed  the  importance  of  finite  grating  depths  for  high 
sensitivity  grating-coupled  waveguide  sensors  by  using  a  rigorous  numerical  analysis.  The  aim  of 
the  present  paper  is  to  analyze  the  influence  of  the  grating  depth  analytically  by  using  a 
perturbational  approach  based  on  the  pioneering  work  of  Kiselev16  in  1975.  Finally,  we  investigate 
the  regions  of  validity  of  the  analytical  model  by  comparing  it  with  numerical  results  obtained 
using  a  boundary  variation  method.17"19 

2.  Wave  propagation  in  a  corrugated  waveguide:  homogeneous  problem 

In  most  grating-coupled  waveguide  sensors  a  coupling-grating  imprinted  into  the  waveguiding 
film  is  illuminated  by  an  external  light  source  while  the  amount  of  light  coupled  into  the 
waveguide  is  being  detected.  However,  as  we  shall  see,  most  of  the  characteristic  sensor  properties 
are  given  solely  by  the  waveguide  geometry  irrespective  of  the  manner  in  which  it  is  being 
illuminated.  In  this  section  we  therefore  start  by  solving  Maxwell's  equations  for  a  waveguide, 
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which  is  left  isolated  without  being  illuminated  from  the  outside.  This  is  referred  to  as  the 
homogeneous  problem. 

A.  Basic  assumptions 

The  system  considered  is  shown  in  Fig.  1.  Here  a  waveguiding  film  with  thickness  dF  and 
refractive  index  nF  is  embedded  between  a  semi-infinite  substrate  with  refractive  index  ns  and  a 
semi-infinite  cover  medium  with  refractive  index  nc .  The  film-cover  interface  is  corrugated  with 
a  sinusoidal  grating  profile  h  of  the  form 

h  (x)  =  <7sin(Kx),  K-2njK,  (1) 

where  a  is  the  grating  amplitude,  K  is  the  grating  wave  number,  and  A  is  the  period  of  the  grating. 

In  the  case  where  no  light  is  incident  upon  the  waveguide,  there  will  only  be  downwards 
propagating  waves  in  the  substrate,  i.e.  waves  with  negative  z  components  of  the  wave  vector,  and 
only  upwards  propagating  waves  in  the  cover.  In  the  film,  however,  both  up-  and  downwards 
propagating  waves  are  present  due  to  reflections  at  the  interfaces. 

We  assume  that  in  each  medium  one  dominating  zeroth-order  wave  with  an  x  component  of  the 
wave  vector  kx  is  present.  Apart  from  that,  diffracted  waves  are  present  due  to  diffraction  off  the 
grating.  However,  the  surface  corrugation  amplitude  a  is  taken  to  be  small  (smallness  to  be 
quantified  later)  so  that  only  the  ±lst  diffraction  orders  in  each  medium  are  taken  into  account. 
This  assumption  is  well  known  from  the  Rayleigh-Fourier  analysis  of  optical  diffraction  in 
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gratings,  which  is  valid  when  the  grating  depth  is  much  smaller  than  the  light  wavelength.  As  a 
result,  a  total  of  twelve  waves  are  included  in  our  analysis,  as  illustrated  in  Fig.  2. 

Like  in  the  conventional  three-layer  waveguide  analysis  without  surface  corrugation,  we  wish  to 
end  up  with  a  mode  equation  that  dictates  which  values  of  kx  (i.e.  which  modes)  are  allowed  in 
the  waveguide.  It  should  be  emphasized  here  that  a  mode  in  our  case  consists  of  a  set  of  twelve 
waves  rather  than  the  traditional  four  waves  considered  for  non-corrugated  waveguides. 

We  restrict  the  analysis  to  cover  only  TE  polarized  waves,  which  is  sufficient  to  illustrate  the  main 
features  of  the  corrugated  waveguide. 

B.  Solution  ansatz 

Since  all  three  media  are  considered  homogeneous,  isotropic,  and  non-magnetic,  we  know  that 
Maxwell’s  equations  have  simple  plane-wave  solutions.  Hence,  for  the  TE  polarized  electric  fields 
we  may  use  the  following  scalar  solution  ansatz  for  the  fields  in  the  substrate,  film,  and  cover 
F 

^ S,F,C  • 


Es  =  'Yjb,  exp[i(kx  +  lK)x-ikzlz\ , 

i=- 1 

1 

*f=Lk  exphX^.  +  lK)x  +  ikFjz\  +  a,  exp|/Ykv  +  IK)x - ikzl  z]},  (2) 

/=- 1 
1 

Ec  =  exp[i(kx  +  lK)x  +  ikzlz] . 

/=- 1 


The  temporal  frequency  term  exp(-icot)  has  been  omitted  for  simplicity,  and  /  is  a  summation 


index  that  runs  through  the  diffraction-order  numbers  -1,0,  and  1.  b,  are  the  amplitudes  of  the 
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three  diffraction  orders  in  the  substrate,  where  the  superscript """  refers  to  downwards  propagating 
waves,  a[  and  a]  are  the  corresponding  amplitudes  of  the  up-  and  downwards  propagating  waves 
in  the  film,  and  If  are  the  upwards  propagating  wave  amplitudes  in  the  cover.  As  mentioned,  kx 
denotes  the  x  component  of  the  zeroth-order  wave  vectors,  which  is  identical  in  all  three  media, 
and  kf  (j  =  S,F,C)  represent  the  z  components  of  the  individual  wave  vectors  defined  by 

Ki  =  ^ln]-(kx+lK)2  .  (3) 

Here  we  have  the  vacuum  wavenumber  k0  =  2x/\  ,  A0  is  the  vacuum  wavelength,  and  nj  are  the 
refractive  indices  of  substrate,  film,  and  cover,  respectively. 

C.  Boundary  conditions 

To  work  out  which  values  of  kx  are  allowed  in  the  waveguide  we  need  to  apply  the  well-known 
boundary  conditions  stating  that  the  tangential  components  of  the  electric  and  magnetic  fields  are 
continuous  across  the  substrate-film  and  film-cover  boundaries. 

(?)  Continuity  of  tangential  electric  field  at  substrate-film  boundary 

By  matching  the  expressions  of  Es  and  Er  from  Eqs.  (2)  at  the  boundary  z  =  —dF  we  obtain  the 
following  three  relations: 


bt  exp  ( ikszldF )  =  a;+  cxp(-ikfdF )  +  a,  exp  [ik'xldF ) , 


(4) 
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where  /  =  -1,0,1.  As  is  seen  here,  the  initial  equation  Es  (z  —  —dF )  —  EF(z  —  —dF )  is  split  up  into 
three  equations  in  Eq.  (4),  which  is  a  result  of  the  orthogonality  of  the  spatial  harmonics 
exp \j(kx  +  IK)  xj  .  Therefore,  the  coefficients  of  each  space-harmonic  can  be  collected  and  set  to 

zero  separately  resulting  in  three  equations  relating  the  nine  amplitudes  b~{  ,  c/;+ ,  and  af  (/  =  - 
1,0,1). 

(ii  )  Continuity  of  tangential  magnetic  field  at  substrate-film  boundary 

From  the  assumption  that  the  electric  field  E  is  TE  polarized,  i.e.  E  =  yE,  and  that  the  media 
involved  are  homogeneous,  isotropic,  and  non-magnetic  it  follows  from  Maxwell’s  equations  that 
the  magnetic  field  H  is  proportional  to  xdzE~zdxE  ,  where  x,  £  are  unit  vectors  along  the  x  and 
z  directions  and  dxz  are  derivatives  with  respect  to  x  and  z,  respectively.  The  tangential  part  of  H 
is  then  found  simply  by  scalar  multiplying  with  the  tangential  surface  vector  t  =x+zdxh  (see 
Fig.  1)  resulting  in 


Tangential  oc  d ZE  ~d xEd xh  ,  (5) 

where  h  represents  the  surface  topography  defined  in  Eq.  (1). 

At  the  substrate-film  boundary  dxh  is  obviously  zero,  which  leaves  us  with  a  simple  matching 
between  d,Es\ _ _d  and  dzEF\  leading  to 


ksz  lb,  exp (ikz  ldF )  =  -k[yf  cxp(-ik)/lr )  +  kz  la,  exp (ikz  ldF ) , 


(6) 
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where  again  /  =  -1,0,1  and  the  coefficients  of  the  spatial  harmonics  exp [i(kx  +  lK)x]  have  been 
collected  and  set  to  zero.  Like  in  Eqs.  (4),  Eqs.  (6)  comprise  three  equations  relating  the  nine 
amplitudes  If ,  c/(+  ,  and  a ~x  (/  =  -1,0,1). 

(Hi)  Continuity  of  tangential  electric  field  at  film-cover  boundary 

By  matching  EF  and  Ec  from  Eqs.  (2),  now  at  the  boundary  z  =  <?  sin  ( Kx) ,  we  obtain: 


I({  u;+  exp  [/&,F;cr  sin  (£*:)[] +  cr;  exp  [-/T^cr  sin ( AA)[]}exp \j(kx  +  IK ) x]  j  = 

(7) 

Z{  If  exp \jkf(J sin  ( Ar)[] exp \f{kx  +  IK )  jc]|  . 

/=-i 

If  we  assume  that  the  grating  is  shallow,  i.e.  a « \/kJzl ,  we  may  use  the  relation 
exp^/k_;;ersin(7^x)J  »1+  ikficr  sin  (Kx)  and  obtain 


1  f  f  kF  <7  1 

y  a,+  <  1  +  — —  [exp  (iKx )  -  exp  (—iKx )]  [  exp  {ilKx)  - 


S 

/=-! 


1 — [exp  (iKx)-  exp  (-/&)]  i  exp  (//£*) 


|l  +  [exp  ( iKx )  -  exp  ( -iKx ) ]|  exp  ( ilKx ) 


(8) 


in  which  the  common  factor  exp (ikxx)  has  been  omitted.  After  expanding  the  sum  in  Eq.  (8), 
constant  terms  (i.e.  terms  that  do  not  depend  on  x)  as  well  as  terms  proportional  to  exp  (-/Ax)  and 
exp  [iKx)  are  collected  and  set  to  zero.  This  results  in  three  equations  relating  the  nine  amplitudes 
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af ,  a,  ,  and  If  (/  =  -1,0,1).  Terms  of  higher  order  being  proportional  to  exp  (-2/ Ax)  and 
exp  (2/ Ax)  are  ommitted,  as  they  relate  to  diffraction  orders  that  are  neglected. 


(iv)  Continuity  of  tangential  magnetic  field  at  film-cover  boundary 
By  using  relation  (5)  to  match  HTangenlial  at  the  film-cover  boundary  we  obtain 


1  f  kF  <7  } 

^  j  1  +  Z'1  [exp  ( /AT )  -  exp  ( -/AT )  ]  >  exp  ( i IKx )  - 


i=-i 


y[  £f,a,  |l  — ^ —  [exp  (/AT)  -  exp  (-/AT)]  j  exp  (//AT)  - 


/=-i 


'^j{kx  +  lK)[a'l  +  at  )  [exp  ( iKx)  +  exp  ( -iKx)  exp  (//AT) 


/=- 1 


k-,<J , 


y[  kc7llf  -j  1  +  — ^ — [exp  (/AT) -exp  (-/AT)]  |  exp  (//AT)  ■ 


/=-i 


A'cr  ^ 


9  ' 

z  /=-i 


y'fkx+lK)  exp  (/AT)  +  exp  [-iKx)  exp(//AT)  , 


(9) 


where  again  the  common  factor  exp(/Aix)  has  been  omitted  and  terms  proportional  to  <r2  have 
been  neglected. 

Like  in  Eq.  (8)  constant  terms  as  well  as  terms  proportional  to  exp  (-iKx)  and  exp  (/AT)  are 
collected  and  set  to  zero  resulting  in  yet  another  three  equations  relating  aF ,  a] ,  and  If  (/  =  - 
1,0,1). 
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D.  Deriving  the  mode  equation 

As  mentioned  above,  each  of  the  four  boundary  conditions  generates  three  equations  which  add  up 
to  twelve  equations  relating  the  twelve  field  amplitudes  cl[  ,  aj ,  Z?,+ ,  and  bj  (/  =  -1,0,1).  However, 
by  combining  Eqs.  (4)  and  (6)  it  is  possible  to  isolate  the  amplitudes  of  the  downwards 
propagating  waves  in  the  form: 


Piai  ■>  ty  y i a j  , 


„  kF,+ks. 

o  _ _ z.l  z.l 

Pi  ~  , 


''z.l 


ks 

Kz,l 


exp (-2 ikFzldF),  y,=  2- 


vj,/ 


'■z.l 


■ks 

Kz.l 


exp 


-i(kF,+ksz,)dF 


(10) 


Using  these  expressions  in  Eqs.  (8)  and  (9)  the  number  of  equations  can  be  reduced  to  six  with  the 
six  unknowns  a[ ,  bF  (/  =  -1,0,1).  After  some  algebraic  manipulation  of  this  system  of  equations 
we  arrive  at  the  following  matrix  equation 


Ae  =  0 , 


(11) 


where  e  =  {a(| ,  <U, ,  ,  b* ,  /y , }  and  the  elements  of  A  are  given  by 
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An=l  +  A>,  A12  =  ~^a-A),  A13=-|<_i(1-^_1),  A14  =  -1, 
A15=|/A,  A16=~kl1,  A^f^a-A),  A22=l  +  A 

A 23  —  A24  —  ~  —  kz0,  A25  —  —1,  A26  —  0, 

A3i=-f  <0(1-A)>  A32=0,  A33=l  +  yff_1,  A34=|^0 
a35  =  0,  a36  =  -1,  A41  =  ^F0(l  -  J30), 

A42=-^[(k^)2+(kx  +  K)K](l+/3l), 

A43  =  §[«-i)2  -  (A  -  *)*]  (1+ /U  a44  = 

A45  =  §  [(C)2  +  (K  +  K)K] ,  a46  =  -f  [(Cl)2  -  (K  -  K)K] , 
A51  =  | [«0)2-^]a+AoX  a52  =  C(1- ax 
a53  =  o,  a54  =  -|[(02-^].  a55=-C  a56  =  0, 

A61  =  -|[(02+^]a+A).  a62  =  o, 

A63  =  Cid-CX  a64  =  |[(C)2+^].  a65  =  o,a66=-Ci. 


To  obtain  nontrivial  solutions  to  Eq.  (11)  the  determinant  of  A  must  equal  zero.  This  results  in  a 
mode  equation  that  can  be  written  in  the  following  form  by  keeping  only  terms  up  to  cr2 : 


D0DlD  1  = 


+TC(C~C) 


(02-(0 


C  \  2 


y 

C(i-A)<,-(i +/AC 


(i+A)-2C 


O.D, 


CT 


--/A 


ADi+ZM,a+A)^ 


Ai  > 


l 


M,  =  - 


2L 


K)  (*5>-^-‘5)+K)  K-^hK) 


T  =  -K,f-2k;X„-(kl,f. 


(13) 
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Here  Dt  (/  =  -1,0,1)  are  the  individual  determinants  of  the  three  groups  of  diffraction  orders  for 


the  non-corrugated  waveguide  ( <7  =  0).  For  example,  D0  represents  the  determinant  of  the  4x4 
matrix  that  would  have  resulted  if  only  the  zeroth  order  waves  ,  a“ ,  /;,} ,  and  b~  had  been 
present  in  a  non-corrugated  waveguide.  Hence,  the  equation  D0  =0  is  equivalent  to  the 
conventional  3-layer  mode  equation  for  a  planar  waveguide  without  grating. 16,20  The  terms 
proportional  to  cr2  on  the  right  hand  side  of  Eq.  (13)  are  small  correction  terms  representing  the 
influence  of  the  grating  on  the  modes. 

E.  Solving  the  mode  equation:  perturbational  approach 

The  mode  equation  (13)  describes  those  values  of  kx  (corresponding  to  propagation  angles  in  the 

waveguide)  that  make  the  original  solution  ansatz,  Eq.  (2),  with  twelve  waves  a  solution  to 
Maxwell’s  equations. 

Usually,  the  solutions  to  the  mode  equation  are  expressed  in  terms  of  normalized  kx  values  called 
the  effective  refractive  index  N  =  k  j k() .  In  general,  a  number  of  distinct  solutions  Nm  exist  for 
each  waveguide  geometry,  with  m  being  the  mode  number.  In  the  following  we  shall  analyze  the 
impact  of  the  surface  corrugation  on  the  solutions  to  the  mode  equation.  Because  cr  is  considered 
small,  we  write  Nm  as  a  sum  of  the  non-corrugated  solution  Nm (a  =  0)  =  /V”'  and  a  perturbation 
term  SNm : 


N_  =  NZ  +  SN,„ 


(14) 
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By  using  Nm  =  N"nc  on  the  right  hand  side  of  Eq.  (13)  it  is  possible  to  recover  a  simplified  mode 


equation: 


n  ~—kF  (kF  -kc  1 

2  '-0  ~’°  Kz-° 


C  \2 


y  K,r-K) 

SUl  -rX/-(l +  '■/)*. 


-d+A)-2  fc 


z,0 


(15) 


To  derive  an  analytical  expression  for  J/Vm  we  expand  Eq.  (15)  to  first  order  in  8 Nm  to  recover 


SN.  =  SNi  +  iSN'i  =  - 


rKo ^ 


2  r 


K +T.K+0 


l=±  1 


vf  exp(-2//c(/r//-)  + 1 
vf  vf  exp(-2  ikFjdp )  - 1 


(16) 


where  <£/Vm  and  <£ZVm  denote  the  real  and  imaginary  parts  of  J/Vm  and  deff  is  the  effective 
thickness  of  the  non-corrugated  waveguide  defined  by2 1 


d  rj  ,  d  p  i 


V  **o  f,o  j 


(17) 


In  most  practical  sensor  applications  light  is  coupled  in  from  air.  In  such  cases  and  when 
M ax { ns ,  nc }  >  (nF  + 1  )/2 ,  which  is  also  mostly  fulfilled,  5 N m  and  8 N m  can  be  written  in  the 
following  analytical  form: 


SNi  = 


f  Ko > 

2  _N2 
,lF  1 V  m 

—2  f 

l  2  J 

k0  N  mdeff 

N: 


+  {ko ) 


kFJ  +  \k 


z.l 


(18) 


12 


SN,„  = 


'  a'' 

X 


<-N2Jn'-r' 

£)■ 

(ko  )(k!,-i)2  +  {k0  )3  (nF  -  rig)  cos2 (k^_xd F 

Nn4eff 

1  1 

7 

to  tT 

+ 

1  _ 1 

2  4 

~(k0)  (nl cos2 (k^dp) 

N=Nm 

(19) 


In  Fig.  3  SN0  and  8N0  are  plotted  versus  the  corrugation  depth  <r  for  the  nanoporous 
waveguide.22'23  It  is  seen  that  SN0  is  negative  implying  that  the  mode  propagation  angle  in  the 
film  8,  defined  by  sin  (6)  =  N0/nF  ,  is  decreased  as  a  result  of  the  surface  corrugation.  As 
expected,  8N0  is  seen  to  be  positive.  The  physical  consequence  of  this  is  that  each  of  the  twelve 
waves  in  the  solution  ansatz  (2)  decay  with  the  damping  factor  exp(-k0J./V0x) .  This  damping  is 

due  to  the  diffractive  out-coupling  of  the  main  wave  into  the  substrate  and  cover  media  through 
the  l  =  - 1  diffraction  orders.  Therefore,  the  radiation  loss  coefficient  of  the  system  is  given  by 

®RAD  ~  k0  ■ 

In  Fig.  4,  8N0  and  SN0  are  plotted  versus  film  thickness  dF  for  a  surface  corrugation  of  <7  =  10 
nm.  Here,  it  is  seen  that  for  large  film  thickness  and  a  film  thickness  close  to  the  cutoff,  SN0  and 
SN0  vanish.  This  is  due  to  the  fact  that  for  a  large  film  thickness,  the  surface  corrugation  has  a 

decreasing  influence  on  the  effective  thickness,  whereas  for  a  film  thickness  close  to  cutoff,  the 
field  penetration  depth  into  the  cover  medium  goes  to  infinity,  so  that  eventually  most  of  the  mode 
power  propagates  outside  the  film  and  therefore  does  not  experience  the  surface  corrugation. 
Mathematically  both  effects  are  described  by  the  1  /  deff  term. 

3.  Coupling  into  the  waveguide:  non-homogeneous  problem 
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In  this  section  we  assume  that  a  monochromatic  plane  wave  is  incident  on  the  waveguide  through 
the  substrate,  see  Fig.  5.  In  order  to  couple  into  the  mode,  we  need  to  launch  the  wave  along  the  - 
1st  diffraction  order,  i.e.  with  an  angle  of  incidence  ds  given  by  ns  sin  (0S )  =  N  -K/k0 ,  so  that 
after  transmission  through  the  substrate-film  interface,  the  incident  wave  gets  diffracted  off  the 
grating  into  the  direction  of  the  a0  wave  components.  The  resulting  solution  ansatz  to  be  inserted 
in  Maxwell's  equations  is  therefore: 


Es  =  c\  exp| i(kx  -  K)x  +  ik^_lz]+'^b^  exp  [i{kx  +  lK)x-ikszlz\ , 

/=- 1 
i 

£,=!{  a[  exp \i(kx  +  lK)x  +  ikzlz]  +  at  exp [i(kx  +  lK)x - ikzJz]} ,  (20) 

/=- 1 
i 

Ec  =  exp [i(kx  +  lK)x  +  ik^z] , 

/=- 1 


where  c\  is  the  amplitude  of  the  incident  field. 


By  following  exactly  the  same  procedures  as  in  Secs.  2C  and  2D,  it  is  straightforward  to  obtain  the 
following  non-homogeneous  matrix  equation: 


A  e  =  c  , 


(21) 


where  A  and  e  are  given  in  Sec.  2D  and  the  driving  vector  c  is  given  by 
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with 


s 

fr,. 
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(23) 


Eq.  (23)  is  readily  solved  by  multiplying  both  sides  with  A  1 : 

e  =  A_1c . 


(24) 


By  keeping  only  the  lowest  order  terms  in  <7 ,  analytical  expressions  for  the  elements  of  e  can  be 
found.  For  example,  the  amplitude  a+{)  assumes  the  form 


a 


o 


<7kz_lk0^nF  nc)  A*’-1c-i 
det(A) 


(25) 


The  total  internal  reflection  of  the  aFl  and  a0  waves  in  the  film  implies  that  |r/(|  |  =  cr0  .  Thus,  the 

l  +l2 

in-coupled  mode  power  is  proportional  to  «0  .  In  Fig.  6  this  in-coupled  mode  power  for  the 

nanoporous  waveguide  versus  N  is  plotted  for  three  different  values  of  surface  corrugation  depth 
<7 .  First  of  all,  it  is  seen  that  for  small  surface  corrugations  a  clear  resonance  close  to  the  non- 


15 


corrugated  effective  mode  refractive  index,  N"c ,  appears.  However,  the  position  of  the  resonance 
peak  shifts  downwards  with  increasing  <7 ,  the  peak  height  decreases  with  increasing  <7 ,  and  the 
half-width  increases  with  <7 . 

These  peak  shape  characteristica  are  largely  described  by  the  l/det(A)  factor  in  Eq.  (25).  This  is 

24  25 

easily  seen  by  Taylor  expanding  det(A)  around  N =  Nm ’ 


det  (A)  =  C,(N-  NJ  +  C2(N  -  Nmf  +  •  •  •,  (26) 

where  Cl  2  are  constants.  Close  to  resonance,  i.e.  N  -7  Nm ,  we  may  keep  only  the  first  order 
term,  after  which  we  get 


det(A)  =  Cx  [ N  -  ( N :  +  SN'J  -  iSNm  ] .  (27) 


i  i— 2 

Thus,  because  the  square  modulus  of  a+  is  proportional  to  det(A)  ,  the  in-coupled  mode  power 
assumes  a  Lorenzian  shape  around  resonance: 


a-(ac+jaC,)]  +(jaC) 


(28) 


where  8Nm  represents  the  shift  in  resonance  and  8Nm  determines  the  peak  height  and  the  half¬ 
width  (=2 8Nm )  achieved  due  to  the  surface  corrugation,  just  as  we  observed  in  Fig.  6. 
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4.  Numerical  verification 


In  order  to  verify  the  validity  of  our  analytical  approach  above,  we  have  analyzed  the  in-coupling 
problem  numerically.  The  numerical  routine  used  is  the  so-called  multiple  interface  boundary 
variation  method.17"19  This  approach  enables  modeling  of  periodic  transmission  consisting  of  an 
arbitrary  number  of  materials  and  interfaces  of  general  shape  subject  to  plane  wave  illumination. 
A  convolution  of  the  exact  waveguide  solution  with  a  multi-layered  boundary  variation  solution 
yields  the  coupling  coefficient  for  light  coupled  into  the  sensor. 

In  Figs.  7a-c  the  results  of  the  numerical  solutions  are  shown  as  dotted  data  points.  As  can  be  seen, 
there  is  good  correlation  between  the  numerical  data  and  the  analytical  curves  for  <7  up  to  ~  80 
nm. 

Regarding  the  peak  shapes,  these  are  further  analyzed  in  Fig.  8  where  8N0  (resonance  position) 

and  8N0  (resonance  height  and  width)  have  been  compared  with  the  analytically  obtained  values. 

It  is  seen  that  the  analytical  model  starts  breaking  down  at  <7  «  80  nm  corresponding  to  about  1/2 
of  the  film  thickness. 

The  8N0  and  8N0  values  may  easily  be  converted  into  corresponding  propagation  angle  values 
89 .  (j  =  S,  F,  C)  in  the  three  media  by  using  that 


89 j  =  8Nq  j tij  , 
89]  =  2 8N'0/nj . 


(29) 
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Here  SO  _  and  80 .  represent  the  shift  in  coupling  angle  and  peak  half-width,  respectively.  Thus, 

for  a  =  80  nm,  corresponding  to  half  the  film  thickness,  the  peak  shifts  and  peak  half-widths 
(measured  in  the  substrate)  are  0.28  deg  and  0.22  deg,  respectively.  Using  the  grating  equation  one 
can  easily  calculate  the  in-coupling  angle,  measured  in  air,  which  is  the  usual  experimentally 
measured  waveguide  sensor  parameter.  For  a  =  10  nm,  a  typical  grating  depth,  the  peak  shift  is 
5.8xlO"3  deg,  which  is  approximately  two  orders  of  magnitude  higher  than  the  resolution  usually 
obtained  in  optical  waveguide  lightmode  spectroscopy.9'26 

5.  Discussion 

Based  on  the  well-known  Rayleigh-Fourier  analysis  we  have  analyzed  the  influence  of  finite 
grating  depths  on  the  TE  mode  properties  in  grating-coupled  optical  planar  waveguide  sensors.  By 
using  a  perturbational  approach,  expressions  for  the  in-coupling  peak  shift  and  peak  width  have 
been  derived.  The  analytical  results  obtained  for  a  specific  waveguide  geometry  were  compared 
with  rigorous  numerical  calculations.  It  was  found  that  for  grating  amplitudes  up  to  approximately 
half  the  waveguide  film  thickness,  the  analytical  model  provides  satisfactory  results. 

In  the  analyzed  example,  it  was  shown  that  a  typical  grating  amplitude  of  10  nm  causes  a  shift  on 
the  order  of  10~3  deg  which  is  two  orders  of  magnitude  larger  than  a  typical  resolution  for  a 
waveguide  sensor  device.  Hence,  the  finite  grating  depth  influence  is  significant  when  performing 
absolute  refractive  index  measurements  or  when  performing  absolute  thickness  detections  of  small 
add  layers,  such  as  proteins,  lipid  bilayers  and  inorganic  film  depositions. 

The  presented  model  can  be  extended  to  cover  TM  mode  propagation  too  and,  moreover,  we 
believe  that  the  method  may  be  used  for  finite  grating  lengths  also. 
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Figure  Captions 

Fig.  1.  Schematic  illustration  of  the  waveguide  structure  considered.  nSFC  are  the  refractive 
indices  of  the  substrate,  film  and  cover  media,  respectively,  a  is  the  grating  amplitude,  t  is  a 
tangential  vector  at  the  film-cover  interface,  A  is  the  grating  period,  and  dF  is  the  film  thickness. 

Fig.  2.  Diffraction  orders  included  in  the  homogeneous  analysis. 

Fig.  3.  Real  and  imaginary  parts  of  the  correction  term  SN0  plotted  against  surface  corrugation 

22  23 

depth  <J  for  a  nanoporous  waveguide.  ~  "  The  parameters  used  are  ns  =  1.22,  np  =  1.57,  «c  =  1-33, 
A0  =  632.8  nm,  A  =  480  nm,  and  dF  =160  nm. 

Fig.  4.  Real  and  imaginary  parts  of  the  correction  term  SN0  plotted  against  film  thickness  dF  for 
a  nanoporous  waveguide.""'23  The  parameters  used  are:  ns  =  1.22,  hf  =  1.57,  nc  =  1.33,  A0  =  632.8 
nm,  A  =  480  nm,  and  <r  =  10  nm. 

Fig.  5.  Wave  vector  scheme  for  the  non-homogeneous  case,  in  which  an  external  wave  with 
amplitude  is  incident  on  the  waveguiding  film  through  the  substrate. 

Fig.  6.  In-coupled  mode  power  versus  effective  refractive  index  N  for  the  nanoporous  waveguide 
(data  given  in  Fig.  3  and  4). 

Fig.  7.  In-coupled  mode  intensity  versus  effective  refractive  index  N  for  the  nanoporous 
waveguide  (data  given  in  Figs.  3  and  4).  Line:  analytical,  points:  numerical. 
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Fig.  8.  SN  and  SN"  versus  grating  amplitude  obtained  analytically  and  numerically. 
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Fig.  5.  Horvath  et  al.,  JOS  A  B 


28 


Fig.  6.  Horvath  et  al.,  JOS  A  B 
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Fig.  7c.  Horvath  et  al.,  JOS  A  B 
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